Oscillations near separatrix for perturbed 
Duffing equation 



O.M. Kiselev, 



Abstract 

A periodic perturbation generates a complicated dynamics close to 
separatrices and saddle points. We construct an asymptotic solution 
which is close to the separatrix for the unperturbed Duffing's oscillator 
over a long time. This solution is defined by a separatrix map. This 
map is obtained for any order of the perturbation parameter. Prop- 
erties of this map show an instability of a motion for the perturbed 
system. 

Introduction 

We consider an equation for the perturbed Duffing's oscillator: 

u" + 2u-2u 3 = ecos(u;t + $ ). (1) 

Here e is a small positive parameter, u and $o are constants. 

The goal of the paper is to construct an asymptotic solution for (pQ) which 
is close to the separatrix for the unperturbed Duffing's oscillator over a long 
time. Formally the words asymptotic solution mean that there exists Eq = 
const > 0, such that a following asymptotic series 

oo 

U(t,e) = y £e n U n (t), (2) 

n=0 

gives a residual as o(s n ), Wn e N, as e G (0, e ) when U(t,e) is being 
substituted into ([I]). 

Our plane is to concentrate on asymptotic solution with a separatrix as 
leading term: 

U Q (t) = tanh(t + t). (3) 
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A similar problem was studied by A. Poincare [T]. He considered a spe- 
cial case of the three body problem. Later V.K. Mel'nikov established that 
the separatrices split under a perturbation and he calculated a gap between 
the separatrices of perturbed equation pj. N.N.Filonenko, R.Z.Sagdeev and 
G.M. Zaslavskii [3j obtained a separatrix map for canonical variables in 
Hamiltonian systems. 

The problem for the perturbation of a separatrix is very important for 
a capture into a resonance. Therefore the problem of dynamics near the 
separatrix was studied in the point of view of the capture. A.I. Neishtadt 
[I] calculated a measure of trajectories which cross the separatrix and are 
captured into a resonance. A.V.Timofeev studied the dynamics near the 
separatrix close to the saddle [SJ. 

A change of the angle variable when the solution crosses the separatrix 
was obtained by J.R.Cary and R.T.Scodjie [6]. D.C.Diminie and R.Haberman 
studied a separatrix crossing near a saddle-center bifurcation and a pitchfork 
bifurcation [7] . Full asymptotic expansions for the problem of the separatrix 
crossing near the saddle-center bifurcation was obtained in [HI E] • 

In this work following new results are obtained. It is constructed the 
map for the dynamics near the separatrix for any power of the perturbation 
parameter. It is shown that the phase shift over a circle near the separatrix 
circle is defined by a term of the order e 2 on the previous circle. The con- 
structed map shows that there exists a manifold of e G (0, e ) for We > 0, 
such that the asymptotic solution has more than A circles near the separatrix 
for VA G N. A Cantor set gives an example of such manifolds. 

This work has following structure. In section [1] the main problem and 
results are presented in formal form. Section [2] contains the derivation of the 
map for the dynamics near the separatrix. In section [3] the consequences of 
the separatrix dynamics are presented. 

1 Main problem and results 
1.1 A trial solution 

Let us consider an asymptotic solution in the form Main problem is to 
construct a bounded asymptotic solution of such type as t G (0, — iVln(e)) 
for any A G N, and e G (0, e ). 
A simplest result is following. 

Theorem 1 There exists two parametric asymptotic solution of (TJP in form 
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{IP and |3j when 

-ln(e)<t<~ln(£). 
Higher-order terms of are 

U n (t) = A-e 2t + B-e- 2t + W~(t), 

as t — > — oo, where 

W~(t) = 0(e- 2nt ) } f--oo, 

and W~(t) has an asymptotic expansion into series of powers of e l ast — > oo, 
such that it does not contains the terms C\e 2t and C^e™ 2 ' where C\ and C2 
are some constants. This means that all parameters for U n (t) are A~ and 
B~ . There fore the parameters of the asymptotic solution are: 



00 



A- = E e n A- 

n=l 

oo 

B- = J2e n B-. 



n=l 



Parameter A~ can be excluded from the solution by time shift in main 
term of the asymptotic solution and $0 hi the perturbation term of ([T]). 
Therefore we take a case A~ = 0. Parameter B~ is a distance between the 
asymptotic solution and the separatrix of unperturbed equation (JTJ). 



1.2 Numeric simulations ; 
studies 

There are two scenarios for a prolon- 
gation of the trial solution. When 
t = 0(— ln(e)/2) the trajectory of 
the trial solution closes to the sad- 
dle point (1,0). This trajectory is 
able to turn to the lower separa- 
trix, which goes to another saddle 
point (—1,0). The different way 
for the trajectory lies near the sep- 
aratrix which goes from (1, 0) to 
(+00, +00). The similar changes 
are possible near the saddle (—1,0). 



challenge for the analytic 
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Figure 1: Trajectories near separatrix on 
the plane (u,u'). 
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Figure 2: The left figure shows the dependence of the oscillating time on e for the 
solution of the Cauchy problem u(0) = 0, u'(0) = 1 for (TjQ) where u = 1, $o = ' K - 
The solutions was studied on the interval t € (0, 100) for the values of the pertur- 
bation parameter e G (0.001,0.1). The right figure shows the scaled structure of 
the peak near e = 0.08. 

These cases are shown on the figure 

m 

Let us concentrate to the trajectories which oscillates between the saddles 
(—1,0) and (1,0). For this we should to find manifolds for the parameters 
of a solution and the perturbation e. These manifolds have a complicated 
structure. For example one can see the dependency of the life-time for os- 
cillating asymptotic solution on the perturbation parameter e. This shown 
on the right figure [3 On left picture [2] one can see the thin structure of the 
peak near e = 0.08. 

This numeric simulations take a challenge for analytic studies. Our goals 
are to calculate the asymptotic solutions, to find a dependency for the tra- 
jectory on the parameters and to give the formulas for the manifolds of the 
parameters for the solutions with an oscillating behaviour. 

1.3 Results 

The asymptotic solution which is denned in Theorem [T] will be used as a trial 
solution and we will prolong it to the large time. The asymptotic behaviour 
of this solution has the same form: 

U n (t)=Ate- 2t + B+e 2t + W+(t), 

as t —>■ oo, where 

W+(t) = 0(e- 2nt ), t^+oc, 
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We will show if Bf < 0, then the solution (j2J), © is bounded when 

-ln(e) < t < -ln(e) 
and when 1 <C t <C — ln(e) the solution has a following asymptotic expansion: 

oo 

u(t, e) = - tanh(0) + £ £ n n n (0), 

n=l 

= t + - ln(e) + - In ( - — £+), 
2 2 v 16 1 ; 

as 6 1 — ► =Foo, where 

u n (9) = a^e ±2e + b^ 2d + w^(9), 

where w^{6) has the same properties as W^if) in Theorem CD 

The following theorem gives conditions for the prolongation of previous 
theorem on iV-circles near the separatrix of the unperturbed equation. 

Theorem 2 The asymptotic solution of (TJP from Theorem[J\may be extended 
on the interval —k ln(e) <C t <C — (fc + ~) ln(e), = 0, 1, . . . , [y] and aas a 
/orm 

oo 

u(t,e) = (-lj^tanhfa) + ^ e n u h n (t k ), 

n=0 

where has the following asymptotic behaviour 

n j 2n p n+1 , 

~ E e2K * fc (E4 E f(«S)t,m < »s(77iarf Jfc )+(i;*)±, im sin(rriwt fc ) 

re=— oo ^ i=0 m=0 

as tfc — > ±oo. The condition of extendability is (— l) m [o"i(m) + Aai(m)] < 
for Vm = 1, . . . , N, N e N. i/ere tae parameter Ao~i(k) is following 

7T 

A<7i(fc) = — 7- — COs(^fc), 

16 cosh(7ra;/2) 

and tfc, and are defined by the recurrent sequence: 

h=t, t k+1 = t k + —ln(e)+ip(k), 



a k (l)=B k , ^(1) = 0, Xk(l) = A k -, 
.o- k+ i(n) + Aa k+1 (n) 



a k (n+l) = -32- 



aAn) + AaAn) 



^{n + 1) = -|(ln (i(ai(n) + A ffl (n))) - ln(2)) + ^(n), 
Xi(ra) = -2, XfcO + 1) = -t^Oi W + A«7i(n))(xfc-i(n) + Axfc-i(n)), 
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Here 



Aa k {n) = {u k n )t >0 , ~ («S)-i,o,o» 
&Xk(n) = (tt*)i 1>0>0 - (O^q.o- 

Corollary 1 T/ie mam term of the asymptotic expansion for u(t , e) depends 
on Bj~, k = 1, . . . , j as tk = 0(1). This means that the asymptotic expansion 
is unstable with respect to a small correction of the parameters. 

Corollary 2 There exist a set of the parameters {A k ~} k *L l , {B^} k x L 1 such 
that the solution has an oscillating behaviour near the separatrices as 5 = 
wln(e) belongs by a Cantor set. 

2 Separatrix dynamics 

In this section theorem [T] is proved in three steps. First step is a construction 
of the asymptotic expansion which is valid near the saddle points u — ±1. 
On the second step we obtain an asymptotic expansion which is valid close 
to the separatrices u = ±tanh(i). On last step we match these asymptotic 
expansions and construct the uniform asymptotic solution, which is valid 
over all domains close to the separatrices and the saddles. 

2.1 Separatrix asymptotic expansion 

Let us to construct an asymptotic solution in the form: 

oo 

u(t,e) = X> n On(*)- ( 4 ) 

71=0 

The leading term of the expansion is a separatrix solution 

Uq = tanh(t) 

of the unperturbed equation: 

US + 2U - 2Ul = 0. 

We will study the asymptotic expansion where the main term is the separatrix 
solution of 

U" n + 2U n - WlU n = f n , (5) 

where 

fx = cos(ut + (j) ), 
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and when n > 1 f n is a polynomial of third order with respect to Uj, Uk and 
Ui for j + k + I = n. 

There are two linear independent solutions of the linearized equation: 



v" + 2v- 6U$v = 0. 



There solutions are: 
1 



cosh 2 (t) ' V2 32 cosh 2 (t) T 4 cosh 2 (t) ^ 8 cosh 2 (t) ' 

Wronskian of the functions v\ and v 2 is equal to unit. 
A general solution of (jHJ) has a following form: 

U n = Vl (t) [ dtf n (t)v 2 (t) - v 2 (t) f dtfnity^t) + A\V\ + B x v 2 . (6) 



sinh(4t) sinh(2t) 



St 



to 



to 



Here t is a constant and Ai, B\ are parameters of the solution. 
Lemma 1 Let 



n-l 



k=—oo 



2(n-l) 



1=0 



■ m=0 



fn~ E e H E E ^l m cos(m^)+i/J im sin(m^) 



t — > ±oo, 



then 



2n p n+1 

E e 2fct (E^ E |^^™t)+&»^M 



i=0 L m=0 



£ — > ±oo. 



To prove this lemma one should use asymptotic behaviours for t>i i2 , substitute 
the formula for the /„ into (jSJ) and integrate. 

Denote = f/^i o,o/4 an d -^n = 16^±i,o,o- These parameters define the 
solution for a large time. The changes of these parameters are: 

AA n = A+-A~, AB n = B+-B~. 

The value of AB\ is defined by Melnikov's integral: 

abi = ^ + - fir = Ho =SS^<ft 

= roo^(cOsM)cOs($ )- 

- sinM) sin($ )) = cos($ ) cosh( ; w2) . 



The following formula gives AA 



i- 



lim 



- sin($ ) ^ y v 2 {t) sxn(ujt)dt - 



2 sin(u;s) — u cos(c<js) 
8{u 2 + A) ' 



,2s 
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2.1.1 Validity of the separatrix asymptotic expansion 

The separatrix asymptotic expansion (j2J) is valid until: 

%±±«1. 

It yields the bounds for the independent variable t: 

1 

e exp(2t) < 1, \t\ < -- ln(e) 

and the parameters of the solution and £>+: 

At < e, Bt < £ VnGiV. 



Lemma 2 There exists an asymptotic solution for equation (QP in the form 
when \t\ < -±ln(e). 

2.2 Saddle asymptotic expansion 

The asymptotic expansion has a form: 

oo 

u(t,e) = ±l + J2e n/2 <(T) r = t + r (7) 

re=l 

near the saddle points u = ±1. 

We use the sigh '+' for the expansion near u = 1 and the sign '-' for the 
expansion near u = — 1. 

The correction terms are defined by the following equation: 

where/i = 0, = cos(u;r — u;r + O ) =t 6(«i = ) 2 and is a polynomial of 
order 3, which is defined by the correction terms with the indexes j, k, I such 
that j + k + 1 = n when n > 3. In a general case is a finite sum of powers 
of e T , sine, cosine and independent variable r. A general formula for n-th 
correction term has a form: 

= a±exp(-2r) +/3±exp(2r) +<(r). (8) 

Here u>J(t) has not the terms Cie 2T and C^e -27 " for VCi,C2 = const as 
r — > ±oo. 
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First and second corrections are: 

uf = af exp(— 2r) + (3f exp(2r); 

uf = af exp(-2r) + ft exp(2r) + i-LL exp(-4r) + exp(4r) - 

1 1 

-n a iPl - IT, a cos ( UJT + 00 - wr ). 

2 4 + 

The n-th correction term is estimated for large r: 

uf = 0(exp(±2n|r|)), n > 1, r — > ±oo. 

The domain of validity for the saddle asymptotic expansions is defined 
by the inequality: 

It yields: 

M < -4 HO- 

Lemma 3 There exists an asymptotic solution in form Ffy as e — » 0, where 
u n (r) has form $\) and |t| <C — |ln(V). 

2.2.1 Matching of the asymptotic expansions 

The matching of asymptotic expansions U(t,e) and u(r,e) yields: 

<*L = °> #£» = °; 
One obtains recurrent formulas for the correction terms with odd indexes: 

-2exp(-2t) = £ 1/2 a+exp(-2r), af = -2, 
T = t + r , r = -ln(e). 

Then: 

±Bf exp(2t) = fte-W exp(2r), (3f = ±B+. 
For higher-order terms as t — > oo one obtains: 

4Afexp(~2t)=e 1 / 2 af n+1 exp(-2r), 

a fn+i = 4^. n ; 

^B+exp(2i) = /3+_i£- 1/2 exp(2T), 
/3 2 Vi = ^+ nGN. 
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The sign of depends on the parameter <ft : 
Pt = ^B+ = S 1 - + icos(0o) 



16 1 1 16 v cosh(7rt,j/2) 

If fii > or the same cos(0o) > — 165f cosh(7ru;/2) then the asymptotic 
solution of (JTJ) goes to infinity, otherwise when cos(0o) < — 16-Bf cosh(7Ro>/2) 
the solution of (JTJ) is close to the separatrix which goes to the left saddle. 

Theorem 3 There exists an asymptotic solution of (TJP, such that this so- 
lution has a form (TJ|) when | ln(e) <C t <C — |ln(e) and /jas /orm (0) u>nen 
-Jln(e) < t < -|ln(e). 

2.3 Lower separatrix branch 

Let us consider the case when < 0, then the main term of an asymptotic 
solution as r — > oo is the separatrix which goes from one saddle (1, 0) on the 
plane (u, u') to another saddle (—1,0). This separatrix is: 

u = — tanh(t). 

The asymptotic expansion, which is close to the lower separatrix has a similar 
form as the asymptotic expansion close to the upper separatrix: 

oo 

u(9,e) = u (9) + Y / e n u n (9). (9) 

n=0 

It is easy to see that the asymptotic solutions (El) has the same form as 
(j4j), where one should use the variable 9 instead of t. 

Equation for n-th correction term has a following form: 

w" + 3u n - Qulu n = f n . (10) 

Here /„ is a polynomial of 3-d order with respect to terms Uj,Uk,Ui, such 
that j + k + I = n, and 

n-l / 2 (n-l) n 

/»~ E e ™[ E ° l E ( F t,i, m cos(muj9)+H^ m sm(muj9) 



k=—oo 1=0 



m=0 



±oo, 



The general formula for the solution of ( {TO}) yields: 



In r n+1 



fc=— oo i=0 



m=0 



E e 2fce E^ E U£, m cos(m^)+^ m sin(m^ 



9 -> ±oo. 



Denote = wi 100 /4 and 6^ = 16tt^ 10 . These parameters define the pa- 
rameters of the solution for the large time. The changes of these parameters 
are: 



Aa n = a n - at, Ab n = b n - 6+ 
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2.3.1 Validity of separatrix expansion 

The lower separatrix expansion is valid until: 

g»n+l j 

It yields: 

sexp(20)<l, |0| < — ln(e) 



Lemma 4 There exists an asymptotic solution of (QP in /orm |0|) as |0| <C 
-|ln(e). 

2.3.2 Matching of lower separatrix asymptotic expansion and asymp- 
totic expansion close to right saddle 

Let fii < 0. In this case one can match the asymptotic expansion in the 
right saddle point and the lower separatrix asymptotic expansion @. 
The matching yields: 

e 1/2 /3 1 + exp(2r) = -2exp(20) and e 1/2 af exp(-2r) = eAaf exp(-20), 

as r — > oo and — > — oo. Here < 0, as a result we obtain: 

2r + - ln(e) + hi(-/3+) = 20 + In 2. 

Or the same 

e = r+iln(e) + iln(- i 9 1 + )-ilB2. 
A substitution of and reductions give a following formula: 

J ln(e) + ln(a+) -2r = ln(e) +2 ln(2) +ln(af ) -2r- - ln(e) -ln(-/?+) -ln(2), 
or 

«+(-/?+) = 8a+ a+ = ictf (-#")■ 

It is easy to see that = —2. A following computation gives formulas for 

(3% and 6^: 

£ 3 / 2 /5 3 + exp(2r) = 4^exp(20). 
lo 

As a result we obtain: 

ln(/#) + 2r = -41n(2) + ln(6+) + 2r + ln(-/5+) - ln(2), 
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or 



Pt = ^bt(-Pt), K = 32- 



The same way leads us to formulas for the higher order terms: 
£ (2n+1)/ V„-iexp(-2r) = 4e»a+exp(-20), 

£ (2n+1)/2 /? 2 + n+1 exp(2r) = l e &+exp(26>), 
(-/?+) " 32°"' ° n " ^ /3+ • 



It is convenient to write out the formulas which connect the parameters 
of the asymptotic expansions near the upper and the lower separatrices: 



a n+i — 



1 

32 



af = ^(Sr + ABx), 
(A+ + AA n )(S 1 - + A J B 1 ), 



+ _ B n+1 + AB n+1 
bn ~ 61 B± + AB 1 ' nGiN ' 

= t + \ He) + \\n ( - + A£x)) - ± ln(2), 



= -- 

y 2 



ln(e) + In ( - l(Sr + AB X )) - ln(2) 



2.4 Neighborhood of left saddle point 

2.4.1 Matching of saddle asymptotic expasion and lower separa- 
trix asymptotic expansion 

The matching of the expansions near left saddle point yields: 

2exp(-2#) = £ 1/2 c^exp(-2a), 

exp(26) = e 1 / 2 ^ exp(2a), 

lo 

e\a{ exp(-26>) = e 3/2 a 3 exp(-2cr). 

As a result one obtain formulas for pthe arameters of the saddle asymptotic 
expansion: 

-20 = ^ln(e)-2<7, 
12 



16 x 

A matching yields formulas for higher-order terms: 

e (2n+i)/2 (a - +iexp( _ 2£r)+ ^- +iexp(2a)) = e "(4a-exp(-20) + ^6^exp(20)). 

There fore one obtains formulas for the coefficients of the asymptotic expan- 
sion: 

Here one can obtain an explicit formula for the changing of coefficient b n . 
For example, 

A6 ' = 6 ' " b ' = ^72) cos ( " >^ + H-#) - H2)) + * 

As a result we obtain the following theorem. 

Theorem 4 If < then there exists an asymptotic solution of (QJ) u>/ien 
|ln(e) <C t — ~ln(e) and has following form |D) when o m ( £ ) ^ * ^ 
-§ ln(e) (0) with sub sign "+" when -\ ln(e) < £ < -| ln(e), (TJj 1 < t < 
— ln(e) and (0) sign w/ien — |ln(e) <C £ <C — lii(e). 

2.5 Next circle near separatrix 

Let us consider an oscillation near the separatrix. Denote parameters of 
asymptotic solution by sequence numbers: 

t = t, 9 = 9, r = r, 
a = a, BZ(0) = B±, 6±(0) = &± 
$(0) = $, 0(0) =0. 

The parameters of the asymptotic expansion for following circles may 
being calculated by formulas: 

if Bf (m) + ABi(m) < and n > 2: 

0(m) = t(m) + ^ ln(e) + \ In (^(^(m) + Ai^m))) - ± ln(2), 

I . . B~^ (m) + AB n +i (m) , , 

6+ m = -32 "+ u . ; "7\ 12 
nK ' Bi{m) + A£?i(m) v ; 



o ' m = 

nv ; 32 



32 

' '^_ 1 (m) + AA»-i(m 



— ( -Bf(m) + ABi(m) ) , 



Bj-(m) + AB^m] 
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if bf(m) + A&i(m) > 0: 

t{m + 1) = 9{m) + - ln(e) + In (— (6+ (m) + A6 1 (m))'| - - ln(2), 
2 \1d / 2 

S- m + 1 = -32 n +^. ^ - w + u ; , 13 
nV ' bj(m) + Afc^m) 

Ar = 3^(^ (m) + A6l(m) )' 
^nO) = -gg a n-i( m ) + Aa n-i(m) + A&i(m)\ 

In particular, 

7r cos($(m)) 



B+(m) = B{(m) + 



16 cosh(^f; 



,-/ x , + / \ ^ cos(0(m)) A . 

t ' w=j;w+ idhfF' (14) 

For given parameter $(0) we obtain following discreet dynamics over a 
near-separatrix circle: 

$(m + l) = u{t{m + 1) - 9(m + 1)) + 0(m) 

= -f (ln(e) + ln(i(6f (m) + A&i(m)) - ln(2)) + 0(w), (15) 

0(m) = u(8(m) — t(m)) + <£>(m) = — 

w(ln(e) + | ln(i(5r(m) + AB^m)) - 
§ln(2))+$(m). 

It is easy to see that <E>(m) and 0(m) have the same order: 

0(m) = -u (- + m) ln(e) + 0(1), 

$(m) = -mu ln(e) + 0(1). 

The shown formulas give a dependency on the sequence {B~}^ =0 for 
dynamics of the solution for times t = 0(iVln(e)) for ViV G N. 

Lemma 5 For \/e G (0, £o)> where e = const > i/iere exists a sequence 
{B~}%L such that Bf(k) < and bf > VJfe. 

The sequence of the lemma and inequalities ffT2l and ffTHj) give 

Theorem 5 For Vfe G (0, Eq), where e = const > there exists the sequence 
{B~}^L such that the asymptotic solution oscillates near the separatrices as 
t = 0(JVln(e)) for all ViV G N. 



14 



3 Discreet dynamical system 



Let us consider the discreet dynamical systems which defines by separatrix 
map (ED, (|13l) . It is convenient to change this system by the following way. 
Let us define 

a k (2n-l) = B k (n-1), a k (2n) = b+{n - 1), nGN; 
Aa k {2n - 1) = AB k (n), Aa k (2n) = Ab k (n - 1); 
Xk{2n-l) = A k (n-l), X k{2n)=a+{n-l), k>2, JiGN; 

Ax fc (2n - 1) = A4(h - 1), A Xk (2n) = Aa k (n - 1); 

ip(2n — 1) = + ncu ln(e), ■0(2n) = 0(n) + (n + ln(e). 

This formulas is prolonged on the next step if (<7i(n)+Ao"i(m))(— l) n > 0. 
The formulas for a discreet dynamical system are: 



Those formulas give a nonlinear discreet dynamical system. The non- 
linearity defines by the terms Aa k (n), Axkin) or the same AB k (n), Ab k (n), 
AA k (n), Aa k {n). 

3.1 Generalized Bernoulli shift 

The separatrix map defines a Bernoulli shift for parameters of this map. 
Formulas (1T21) and (|T3l) show that the a k {n + 1) depends on a k+ 2(n). It 
means that (n)-th correction on the (m + l)-th oscillation depends on (n + 2)- 
th correction on the m-th oscillation near the separatrix. This is typical 
property of the Bernoulli shift. This property shows the loss of the accuracy 
for the approximation of the motion and the instability the solution with 
respect to initial data. 



cr k (l)=B k , X k(l) = A k , V(l) = $, 
if (<7i(m) + A(x 1 (m))(-l) m > 0, Vm < n, then 




1>(n + l) = ~(hx (^(<n (n) + Am(n))) - ln(2)) + ^(n), 

Xi(2n) = 2, X i(2n-l) = -2, 
1) = -— (xA-i(n) + Ax*-i(n) )(«7i(n) + A<ri(n)). 
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3.2 Cantor manifold 

When we study the asymptotic solution we say that there exists So > 
such that for some e G (0,£o) there exists studying asymptotic solution. In 
this subsection we turn the study by the other side. We concentrate on 
the problem for the structure of the manifold e G (0,£o) which gives the 
oscillations near the separatrices. 
The following condition 

cos(ij)(n) — 10 — - — ln(e)) < — 16ci(n) cosh (16) 

defines the possibility to prolong of the discreet dynamical system on (n+1)- 
th step. Define S = — ula(e), then S G (So, oo) where So = — u\n(eo) and 
< So < 1. Let <Jx(n) is such that 

arccos(— 16cri(n) cosh(— )) — ip(n) = ±— + 2nk. 

Then (1161) defines the parameter S 

r uo(n — 2) , . . , 
Vk> -4 J -]n(e )] + l. 



Let 



then when 



- S G (-7T + 27ik, 7T + 2nk) G (S , oo), 



-5 < h and 5 > - + 2irk 



2 3 ' 2 3 

the discrete dynamical system is prolonged by the next step. 
Let us consider sequences a G N such that 

— 16cr(Z) cosh(— ) > 1, as I > n a , I < 3n a — 2 

and n a+1 = 3n a — 2, 

arccos(-16ai(n Q+1 ) cosh(^)) - ip(n a +i) = ±tt + Znk, 

Mk > [ U{na+ A l ~ l) ln(e )] + 1. 

Let e G (0, £o), then the set of S, for which the discrete dynamical system is 
prolonged, is the Cantor set on any (— n + 2n, kn + 2nk). 
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